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ABSTRACT 

The nucleosynthesis and other observable consequences of a nova outburst 
depend sensitively on the details of the thermonuclear runaway which initiates 
the outburst. One important source of uncertainty in our current models is the 
nuclear reaction data used as input for the evolutionary calculations. We present 
preliminary results of the first analyses of the impact on nova nucleosynthesis of 
all reaction rate uncertainties considered simultaneously. 



1. Nuclear Uncertainties in Novae 

Observations of nova outbursts have revealed an elemental composition that differs 
markedly from solar. Theoretical studies indicate that these differences are caused by the 
combination of convection with explosive hydrogen burning which results in a unique nu- 
cleosynthesis that is rich in odd-numbered nuclei such as 13 C, 15 N, and 17 0. Such nuclei 
are difficult to form in other astrophysical hosts. Many of the proton-rich nuclei produced 
in nova outbursts are radioactive, offering the possibility of direct observation with 7-ray 
instruments. Potentially important 7-ray sources include 26 Al, 22 Na, 7 Be, and 18 F. 

The observable consequences of a nova outburst depend sensitively on the details of the 
thermonuclear runaway which initiates the outburst. One important source of uncertainty 
is the nuclear reaction data used as input for the evolutionary calculations (Starrfield et al. 
1998). A number of features conspire to magnify the effects of nuclear uncertainties on nova 
nucleosynthesis. First, the similarity of the nuclear burning and convective timescales in 
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novae results in nuclear burning which is far from the steady state which typifies quiescent 
burning. Further, many reactions of relevance to novae involve unstable proton-rich nuclei, 
making experimental rate determinations difficult. Finally, for hydrodynamic conditions 
typical of novae, many rates depend critically on the properties of a few individual resonances, 
resulting in wide variation between different rate determinations. As a result, statistical 
model (Hauser-Feshbach) calculations, which are employed with great success in other venues 
for a large number of reactions (Rauscher et al. 1997), are unreliable for many rates under 
nova conditions. As we will demonstrate, these factors greatly restrict our ability to make 
precise predictions of the nucleosynthesis of novae. 



Though analysis of the impact of variations in the rates of a few individual reactions 
has recently been performed using one-dimensional hydrodynamic models (Jose et al. 1999), 
analysis of the impact of the complete set of possible reaction rate variations in such hydro- 
dynamic models remains computationally prohibitive. We therefore begin by examining in 
detail the nucleosynthesis of individual zones, using hydrodynamic trajectories (temperature 
and density as a function of time) drawn from nova outburst models. Such post-processing 
nucleosynthesis simulations have been a common means of estimating nova nucleosynthesis 
(see e.g., Hillebrandt & Thielemann 1982; Wiescher et al. 1986). For this presentation we 
are using a hydrodynamic trajectory for an inner zone of a 1.25M Q ONeMg WD, similar to 
that described by Politano et al. (1995). These calculations were performed using a nuclear 
network with 87 species, composed of elements from n and H to S, including all isotopes 
between the proton drip line and the most massive stable isotope. 

To investigate the extent to which nuclear reaction uncertainties translate into abun- 
dance differences, we use a Monte Carlo technique which simultaneously assigns a random 
enhancement factor to each reaction rate in the nuclear network. Monte Carlo methods have 
been employed with success in the analysis of Big Bang nucleosynthesis (Smith et al. 1993), 
but this is their first application to other thermonuclear burning environments. This is also 
the first analysis of nova nucleosythesis to examine the impact of all the reaction uncer- 
tainties simultaneously. The reaction rate enhancement factors are distributed according to 
the log-normal distribution, which is the correct uncertainty distribution for quantities like 
reaction rates which are manifestly positive (Smith 1991). The probability density function 
is given by 



where a and (3 are the (logarithmic) mean and standard deviation. Unlike the normal 
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Monte Carlo Estimates of Uncertainties 






Fig. 1. — Histograms showing relative deviations in abundance. Y Center, Y Upper, and Y 
Lower correspond, respectively, to the central abundance and the upper and lower error bars 
from Fig. 2. The dotted curves are normal distributions with the same (arithmetic) mean 
and standard deviation as the Monte Carlo distribution. The dashed curves are log-normal 
distributions with (logarithmic) mean and standard deviation from the Monte Carlo. 



(Gaussian) distribution, where confidence intervals are symmetric in the additive sense, the 
confidence intervals for the log-normal distribution are symmetric in a multiplicative sense. 
For small uncertainties (< 20%), the difference between the log-normal distribution and 
the normal distribution is small. However, for uncertainties of larger sizes, such as those 
encountered in this problem, the difference is important. For example, assuming a normal 
distribution with a relative uncertainty of 50% implies that there is a 5% chance of a negative 
reaction rate. 

The choice of uncertainty for each reaction rate is crucial for this analysis. For /3-decays, 
we have used the experimental uncertainties (Tuli 2000). For many of the other rates im- 
portant for nova nucleosynthesis, the uncertainties in the reaction rate are unknown or are 
likely to be dominated by systematic effects such as missing resonances. We have therefore 
assigned uncertainties by category. Rates calculated by Hauser-Feshbach methods as well 
as rates whose measurement require radioactive ion beams were assigned uncertainties of 
~ 50% (P = ln(1.5)). For rates measured with stable beams, we have assigned (3 = ln(1.2). 
Note that relatively few reactions, especially among unstable nuclei, have measurement un- 
certainties this small. We have chosen these small uncertainties to obtain a conservative 
estimate of the impact of nuclear uncertainties on nucleosynthesis predictions. 

The histograms drawn in Figure 1 show the abundance distributions which result from 
the Monte Carlo analysis (992 trials) for two representative nuclei, 13 C and 22 Na. From 
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Fig. 2. — Final Abundances (elapsed time = 4.7 x 10 5 sec after peak). 

distributions such as these, mean values and 90% confidence limits are extracted. The 
normal and log-normal distributions with widths drawn from these Monte Carlo abundance 
distributions are also plotted in Fig. 1. Note the similarity between the normal and log- 
normal distributions for 13 C, which has a standard deviation of 12%, while these distributions 
differ markedly for 22 Na, which has a standard deviation of 32%. 

Figure 2 shows the mass fractions of each species at the end of the simulation, 4.7x 10 5 sec 
after peak temperature. Because of the long time which has elapsed, the unstable proton-rich 
nuclei have decayed, reducing their abundances to less than 10~ 20 . The error bars displayed 
in Fig. 2 are the 90% confidence intervals which result from the Monte Carlo analysis. As 
evidenced by the error bars in Fig. 2, the impact of even our conservatively chosen variations 
in reaction rates on the nucleosynthesis is large. While broader conclusions will require 
similar analysis of an entire model, a number of interesting points can be drawn from the 
analysis of this single hydrodynamic trajectory. First, for reaction rate variations of the size 
considered here, the impact on the rate of energy production is small. At the 90% confidence 
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level, variations in the amount of hydrogen consumed and helium produced represent ~ 10% 
variation in the thermonuclear energy released. This is essential if post-processing techniques 
are to prove useful. Second, for the most abundant CNO isotopes, such as 13 C and 15 N, 2a 
variations of 5-20% are common, though some less abundant CNO isotopes show larger 
variations. For example, the 90% confidence interval for the mass fraction of 17 has a lower 
limit of 2.7 x 10~ 5 and an upper limit of 5.3 x 10 , a span of a factor of 20. Third, for many 
of the nuclei in the NeNa and MgAl cycles, including the 7-ray source nuclei 22 Na and 26 Al, 
the 90% confidence interval spans variations of a factor of 3. This represents nearly a factor 
of 2 difference in the distance from which novae may be observable with 7-ray telescopes. 

3. Determining the impact of individual reactions on the nucleosynthesis of 

individual isotopes 

While the use of the Monte Carlo described in the previous section allows us to estimate 
the uncertainty in our nucleosynthetic calculations caused by nuclear reaction uncertainties, 
it does not indicate which reactions influence the production of an individual isotope. In 
order to determine this connection, we re-examine our Monte Carlo trials to find correlations 
between the abundance variations and the variations in all of the included reactions. Given 
the wide range of variations, we test for correlations in the logarithm of the abundance change 
with respect to the logarithm of the reaction variations. Figure 3 displays some representative 
correlations between reactions rates and the abundance of 22 Na. Since we simultaneously vary 
all the included reactions rates, the resultant plots show some scatter, however statistically 
significant correlations are still achieved for a number of reactions. Fig. 3a, Fig. 3b, and 
Fig. 3d show that increases in 21 Na(p, 7) 22 Mg, 22 NaQo, 7) 23 Mg, and 23 Na(p, 7) 24 Mg resulted in 
reduced abundances of 22 Na, while Fig. 3c and Fig. 3e show that increases in the reaction rate 
of 20 Ne(p, 7) 21 Na and 23 Na(p, a) 20 Ne increased the abundance of 22 Na. Interestingly, Fig. 3f 
demonstrates that, for reaction rate variations of the size considered here, the abundance 
of 22 Na is essentially independent of the rate of 21 Ne(p, 7) 22 Na, even though this reaction 
in the primary production channel for 22 Na. Much larger variations in the reaction rate of 
21 Ne(p, 7) 22 Na would necessitate reexamination of this conclusion. 

Jose et al. (1999) have recently performed a set of spherically symmetric nova outburst 
simulations in which they vary individual reaction rates in order to examine the impact of 
these rates on the production of 22 Na and 26 Al. They found that variations in the rates of 
21 Na(p, 7) 22 Mg and 22 Na(j», 7) 23 Mg produced significant variations in the 22 Na abundance, 
with 23 Na(p, 7) 24 Mg having a smaller effect. Our results confirm these conclusions, however 
we find that the reaction 20 Ne(p, 7) 21 Na has effects similar in strength to the former two and 
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Fig. 3. — Correlations of variations of individual reactions with the abundances of 22 Na The 
three lines portray the most probable values of the correlation bounded by its 90% confidence 
limits. 
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that 23 Na(p, a) 20 Ne has an effect similar in strength to the latter. Examination of these rates 
within hydrodynamic simulations of novae would present an excellent test of the predictive 
power of the Monte Carlo analysis. 

4. Conclusions 

We have performed analysis of the effects of reaction rate uncertainties on nova nucle- 
osynthesis. Our Monte Carlo analysis is the first such analysis that simultaneously examines 
the impact of all of the uncertainties. The large uncertainties in nova nucleosynthesis that 
we find constrain our ability to make detailed comparisons between theoretical models for 
the nova outburst and astrophysical observations to a degree which is often ignored. Fur- 
ther, knowledge gleaned from this analysis of which reaction rate uncertainty dominates the 
uncertainty in a given isotope helps prioritize new experimental measurements. Only with 
improved knowledge of these uncertain rates, both experimental and theoretical, can we 
provide tight constraints on the nova outburst from its nucleosynthetic products. 
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